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Abstract 

We present a combined analysis of the space- and time-like electromagnetic pion form factors 
in light-cone perturbative QCD with transverse momentum dependence and Sudakov suppres- 
sion. Including the non-perturbative "soft' QCD and power suppressed twist-3 corrections to the 
standard twist-2 perturbative QCD result, the experimental pion data available at moderate ener- 
gies/momentum transfers can be explained reasonably well. This may help towards resolving the 
bulk of the existing discrepancy between the space- and time-like experimental data. 
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Introduction. The electromagnetic (e.m.) form factor of hadrons are important physical 
observables that play a key role in understanding the transition from the perturbative to 
the non-perturbative behavior in particle physics. The space- and time-like e.m. pion form 
factors F n and G n , respectively, are specified through the following matrix elements: 

e(P' + P), F^Q 2 ) = (tt^P') I JJT(0) | ^(P)) , 

e(P>-P),G n (Q 2 ) = (n+(P>)n-(P)\j-(0)\0) , (1) 

where J^ m is the e.m. current, and P = (Q/y/2, 0, T ) and P' = (0, T ) are, respec- 

tively, the initial and final state external light-cone pion 4-momenta in the Breit-frame. For 
the space-like momentum transfers, q 2 = (P' — P) 2 = —Q 2 < 0, whereas for the time-like 
momentum transfers q 2 = (P' + P) 2 = Q 2 > 0. 

Theoretical predictions based on standard "asymptotic" QCD rely on collinear factoriza- 
tion [1] that lead to the celebrated quark counting rule, {F, G} 7T (Q 2 ) ~ l/Q 2 [2'j]. Naively, 
one may then conclude that at high enough energies/momentum transfers the space- and 
time-like form factors are essentially of the same magnitude. However, at the present exper- 
imentally accessible energies the reported pion form factor results differ significantly, with 
the time-like results js-5] being up to a factor of four more than the space-like results js]. 
Efforts to explain the above difference with conventional Vector Meson Dominance (VMD) 
6) and perturbative QCD (pQCD) lead to the general conclusion that the time-like and the 
space-like data are inconsistent with each other. The purpose of this letter, is to show a 
possible "clean" scenario where the above difference could be reconciled with the standard 
treatment of the parton transverse momentum dependence (TMD), the sub-leading twist-3 
contributions and the so-called "soft" QCD corrections. However, in dealing with the par- 
ton picture in pQCD, one should naturally be aware of the fact that in reality there are 
additional difficulties with hadronization and other final state interactions, and resonances. 
Our approach assumes that these effects for Q 2 ^> Aq CD are comparatively small (as data 
suggest) and would not come in conflict with our predictions that account for the largeness 
of the existing discrepancy between the space-and time-like pion form factor data. 

The charged pion form factor can be written as [3] {F, G} 1T {Q 2 ) = {F, G} s ° ft (Q 2 ) + 
{F, G}^ vd (Q 2 ). The factorizable hard part {P, G}^ vd (Q 2 ) is calculated using light-cone 
pQCD with explicit TMD of the constituent valence partons; whereas, the non-factorizable 
soft part {P, G} s ° ft (Q 2 ) is modeled using QCD sum rules (QCDSR) via local quark-hadron 



duality. 

Parametrically, both the soft and higher twist contributions to the form factor are ex- 
pected to be small at large momentum transfers compared to the leading hard (twist-2) 
contributions due to the relative 1/Q 2n suppression. Despite this, their contributions turn 
out to be unnaturally large at moderate range of energies. In this paper, for the first time, 
we show that the twist-3 corrections to the time-like pion form factor are very large and es- 
sentially account for the bulk of the observed discrepancy between theory and experimental 
data. Note that the first attempt to explain both the space- and time-like data in the context 
of pQCD includes only the twist-2 effects [8|. However, the present consensus is that the 
twist-2 effects are much too small to ex plain t 
must use appropriate Sudakov factors 



re form factor data [71. l9Hl2j . Furthermore, one 



13-17] to suppress the kinematic enhancements that 

nnr 

may invalidate factorization. The advantage of such a modified 'fc-r'-factorization |13l. Il4j. |18| 
approach is the elimination of large logarithms in the hard kernel through the TMD of the 
valence partons. This extends the range of applicability of pQCD down to very moder- 



ate range of energies and has been wide 
especially, to exclusive S-meson decays 



v applied to inclusive and exclusive processes, and 



applied to 1 
l4U,ll9|. 



Factorized pQCD. We now present the essentials of our calculations. The dominant contribu- 
tions come only from the leading order (LO) Fock state, i.e., a qq valence quark configuration 
with one hard gluon exchange in the scattering kernel sandwiched between 2-particle wave- 
functions/distribution amplitudes (DAs). One of four diagrams contributing to each of the 
Born amplitudes tc-/* — > tt and 7* — > 7r + 7r~ is shown in Fig. [TJ The other diagrams corre- 
spond to allowing the gluon to interact on the other side of the photon vertex and allowing 
the photon also to couple to the other valence quark. The higher Fock state contributions are 
neglected being suppressed by higher powers of 1/Q 2 - Since, we are only concerned with the 
leading Fock states in the scattering kernel of the pion, we must consider only the 2-particle 
pion DAs for our analysis neglecting the multi-particle components. Nevertheless, one can 
show that for the 2-particle twist-3 DAs are not independent of the 3-particle twist-3 DA, 
being related by QCD equations of motions. To next-to-leading order in conformal twist 
there is just one 2-particle twist-2 collinear DA </> 2;7r (a;, //) with an axial-vector structure, and 
two 2-particle twist-3 collinear DAs, one with a pseudo-scalar structure ^^(x,^) and the 
other with a pseudo-tensor structure ^^(x, fi). They can be derived from light-cone QCD 
sum rules (LCSR) and are usually expressed as truncated conformal series expansion over 
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FIG. 1: LO representative matrix elements contributing to (a) the space-like F 7T (Q 2 ), and (b) the 
time-like G n (Q 2 ) pion form factors. The blobs represent the pion wavefunctions PV 



Gagenbauer polynomials 20H22j . Their asymptotic forms are given by 



i ( as ) 
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x(l-x), (2) 



where /„- ~ 131 MeV is the pion decay constant and x is the longitudinal momentum fraction 
of the valence partons. The intrinsic TMD of the total pion wavefunctions is modeled via 
the Brodsky-Huang-Lepage (BHL) prescription [23] having the general impact "6"-space 
representation: 



V t -n(x, b, //, M q ) = A t . n i;7r (x, /i) exp 



exp 



b 2 x(l 



x 



t = 2,3 



(3) 



x(l — x) 

where 0t ;7r (x,/i) is one of the twist-2 or twist-3 non-asymptotic collinear DAs at any given 
scale fi. Such DAs satisfy the Efremov-Radyushkin-Brodsky-Lepage (ER-BL) evolution 
equation e.g., the twist-2 DA is given at the LO by the following non-asymptotic ex- 
pression, in terms of Gagenbauer polynomials Cn 2 (2x — 1): 



-4 7 i 0) /9 



+ 0(a s ) 



(4) 



n=0,2,4,- 

where a s is the standard (two loop) MS QCD coupling with Aq CD = 0.2 GeV, a^'s are 
the moments of the DA, depicting the genuine non-perturbative inputs, and jn^s are the 
corresponding standard LO anomalous dimensions. A compilation of the numerical values of 
the Gagenbauer moments as well as the LO RGE behavior of the various non-perturbative 



parameters of the twist-2 and twist-3 DAs can be found in 12j, [22|, normalized to the mass 
scale [Aq — 1 GeV. 



The BHL Gaussian parameters A t -^ and /3t- n are fixed using phenomenological constraints 
from 7T° — > 77 and n — > fiu^ decays (see, e.g., [l2|), and M. q ~ 0.33 GeV is the constituent 
(q = u, d) quark mass, introduced to parameterize the QCD vacuum effects. These pa- 
rameters could be additionally constrained using a combined analysis of data from lattice 



simulations and from experiments like C 
production. However, since the analysis 



EO, BaBar and FermiLab E791 diffractive dijet 



12] showed that the sensitivity to the model DA 



parameters is less than 5%, while the experimental error bars are much larger, we refrain 
from doing such a involved analysis at the moment. With the availability of higher quality 
data in future such a systematic combined analysis may provide important constraints to 
our results. 

Next, using the TMD modified factorization ansatz in the operator convolution form 
{F, G}^ ard ~ V n <8> .Mlo <8> "^V, where the LO matrix elements A^lo are diagrammatically 
represented in the FigJT] and £§> represents the phase space integration, one can obtain the 
pQCD contribution to the hard form factor in a standard way up to twist-3 corrections, 
given by {F, G}^ ard (Q 2 ) = 5{F, G}[ twist2) (Q 2 ) + 5{F, G}[ twist3) (Q 2 ) where, 



5{F,G}^ wist2 \Q 2 ) 



64tt 



Q / dxdy 



3 ./a jo 



bidb 1 b 2 db 2 oc a (t) [±x V 2 ^{x, b\) Vz^(y, b 2 )] 
x H±{x, V, Q, h, b 2 ) S t (x) exp [S(x, y, h, b 2 , Q)} ; (5) 



5{F,G}l twist3) (g 



^— ^ nl i dxdy j bidb 1 b 2 db 2 a s (i) 
3 Jo Jo 



x 



xVUxM)VUyM) 



d.V&ix, 6 X ) Vi in (y, b 2 ) + lp^( Xi h) V&(y, h) 



x r H±(x,y,Q,bi,b 2 ) S t (x)exp [S(x,y,b 1 ,b 2 ,Q)] 



(6) 



In the above equations, "+" and "— " correspond to the space-like and time-like cases, 
respectively, V t . )7r (x,b) = P t;7r (x,6, l/b,A4 ud ) and t = max(y / xQ, 1/bi, l/b 2 ) is related to the 
factorization scale. The so-called "chiral" parameter /j, n arises from the standard definitions 
of the twist-3 collinear DAs defined at a suitable low energy scale, fi n (fio ~ 1 GeV) = 
m 2 /(m u + nid) ~ 1.7 GeV 2^|. However, in the context of intermediate energies, fi w is 
usually taken to be slightly lower f» 1.3 — 1.5 GeV which is consistent with fits to the B — > ic 



transition form factors 



fly, 



171 . Il9j , %PT estimates 2l|, |24j and the moment calculation 



applying QCDSR [25|. Since, the twist-3 results can be somewhat sensitive to this parameter, 
here we use /jl^ — 1.5 ± 0.2 GeV and indeed show that its variation contributes to a large 



uncertainty in the time-like region. The hard kernels ~H± could be expressed in terms of the 
standard Bessel functions K (9), I (6), H^(6) = J Q (9) + iY (9) and J o (0): 



U + (x,y,Q,b 1 ,b 2 ) 



K (y/xyQb 2 ) 

x [e(b 1 -b 2 )K (^Qb 1 )I (V^Qb 2 ) 
+ 6(b 2 - h^oiVZQhMVZQbi)] 



(7) 



H-{x,y,Q,bi,b 2 ) 



xy Qb 2 ) 

0(&i - b 2 )H^ 1 \^Qb 1 )J (V^Qb 2 
+ 6(b 2 - b 1 )Hg ) (V^Qb2)MV^Qbi) 



x 



(8) 



The Sudakov factor S(Q) and the jet function S t (x) are introduced to organize to all orders 
the large double logarithms a s \n 2 kT {hr is the generic transverse parton momenta) and 
a s ln 2 x, respectively, that arise from radiative gluon effects and may otherwise invalidate 
perturbative factorization. Such resummations result in the natural suppression of possible 
non-perturbative and kinematic endpoint enhancements of the scattering kernel, thereby, 
improving convergence and making perturbative evaluation self-consistent. For their explicit 
expressions, one is referred to 12Hl7|. A few comments regarding our factorized results 
(Eqs. [5] and [6]) are now in order: 

1) Here, we have presented a LO analysis of the hard kernel which is apparently gauge 
dependent (light-cone gauge), arising from the contribution of the single gluon propagator. 
However, in 26] it was shown that for the tc-/* — > 7 transition form factor, the gauge 
invariance of the hard kernel is a consequence of the gauge- dependence cancellation between 
the quark level diagrams of the full QCD and effective diagrams of the pion wavefunction, 
order by order in perturbation theory using the principle of mathematical induction. In this 
way, the hard kernel and the resulting predictions from the /c^-factorization turn out to be 
gauge-invariant to all orders. The above reference also claims that such an approach could 
be extended to other elastic and transition form factors, at least up to the level of NLO 
corrections. 

2) Our result for the hard form factor depends on the renormalization/factorization scale 
which is typical of all fixed order calculations. The Sudakov factor that resums a certain 
class of radiative soft-gluon contributions to all orders in perturbation theory is inherently 



factorization scale dependent, while the LO hard kernel that is used to evaluate the hard 
form factor depends on the renormalization scale through the running of a s . In this case, 
the scale dependence is minimized by adhering to a fixed prescription with the renormal- 



ization/factorization scale set to the momentum transfer Q 26|, |27|]. It is, however, believed 
that a systematic higher order calculation can eventually absorb this scale dependence. 

3) There may be a simple rationale why the TMD factorization is expected to work at the 
level of 1/Q 2 power suppressed corrections, although a more rigorous proof is beyond the 
scope of this paper. Firstly, note that the "active" soft gluons which may arise e.g., from the 
3-particle twist-3 DA that probe the hard kernel, bring about additional power corrections. 
Compared to the 2-particle twist-3 corrections considered in this work, the 3-particle twist-3 
corrections is not chirally enhanced (there is a large parametric enhancement from \x v in the 
definition of the 2-particle twist-3 DAs, which brings about a sensitivity to the chiral scale), 
and should be numerically small. Secondly, the rest of the "long-distance" soft gluons that do 
not interfere with the hard kernel may break the TMD factorization. However, in the large 
Q 2 limit, a hadron tends to have a small "color-dipole" due to the Lorentz contraction and 
the Sudakov suppression. Such gluons can not probe the small "color-dipole" configurations 
of qq within the hadronic bound state, and their effects cancel each other. This is the so- 
called "color transperancy hypothesis" . With this assumption, one only needs to care about 



collinear gluon effects and their factorization. Using similar arguments, the authors in [28 1 
have explicitly proven TMD factorization at the twist- 2 level and collinear factorization at 
the twist-3 level. Hence, it is our assumption that the approach presented in the above 
reference can even be straightforwardly extended to include the twist-3 TMD factorization. 



Soft QCD. Next, following [7J, we include the soft 
Local Duality (LD) for the space-like form factor [2 



eynman mechanism) contribution via 



^° ft W 2 )lLD - 1 - (1 V4s %?)3/2 ' ^ 

where so ~ 0.68 GeV 2 is the duality interval for higher excited and continuum thresholds 
which is very naturally almost the "middle" between pion mass m 2 m and that of the 
A\ resonance m\ ~ 1.6 GeV 2 . The VMD models and %PT predictions are not expected 
to work beyond ~ 1 GeV, while standard pQCD with only twist-2 operators completely 
fails to explain the available experimental data. The soft contribution, on the other hand, 



is significantly large at moderate energies 



7| and so are the twist-3 power corrections [12] . 
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FIG. 2: The relative magnitudes of the soft (double-dot black lines), twist- 2 (thin solid blue lines) 
and twist-3 (thick solid red lines) corrections to the pion form factor. The twist- 2 and twist-3 
corrections without including the pre- factors are also displayed. 

However, both the soft and the twist-3 corrections are expected to fall off rapidly as ~ 1/Q 4 
for large Q, so that asymptotically (Q — > oo) one recovers the rigorous leading twist-2 
contributions ~ 1/Q 2 which dominate the form factor. This aspects is clearly revealed 
through our analysis (see, Figs. [2] and E}. 

To extend the analysis to the time-like region, one may analytically continue Eq. M from 

n 

the space-like region. Using such a model ansatz, the authors of [30| were able to show for 
the first time a much larger contribution to the form factor in the time-like region than in 
the space-like, and hence were partly able to resolve the bulk of the discrepancy for large 
Q. However, this gives rise to a single pole at 4s ~2.71 GeV 2 which does not correspond to 
any of the real physical bound states or resonances (e.g, p,u,...) seen in the time-like data. 
In fact, the observed spectrum around 2.7 GeV already appears to be rather "smooth" 
and well above the resonance region (below irui threshold). Hence, rather than trying to 
reproduce the actual time-like data, including the various bound states and resonances, we 
try to explain the continuum contribution with a smooth G s ° ft which has the same leading 
1/Q 4 dependence as F* oft under analytic continuation. Thus, we choose the same form of 
the time-like soft factor, i.e., G s f\Q 2 ) = F* oft (Q 2 ) + 0(1/ Q 6 ) for large Q. Here again 
we should stress that our analysis is entirely based on the assumption that the physically 
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observed low-lying resonances would not spoil the continuum contribution which appear as 
"superposed peaks" on a continuum spectrum. 

Results and Discussion. To this end, it may be notable that the rather ad hoc incorporation 
of the soft part from QCDSR have no a priori correspondence with the hard parts, and 
therefore, may lead to the possibility of some double-counting between the respective soft 
and hard contributions at the intermediate regime. However, such double counting could 
partly be removed by imposing the the vector Ward-identity {F, Gj-^O) = 1. Following the 
argument detailed in [3, Q] , we introduce appropriate power correcting pre-factors to restore 
the Ward-identity, and hence we arrive at our final expression for the space- and time-like 
pion form factors given by 



The above pre- factors of 5{F, G}r ' and 5{F, G})T ' ensure a "smooth" matching of 
the different power-law Q 2 behavior between the soft and the hard parts that preserve the 
gauge invariance condition {F, G}^ ard (0) = 0. In principle, this vector Ward identity can also 
be achieved with larger n values in the Q 2n / ((2so) n + Q 2n ) factors in front of the hard part. 
However, as n — > oo, the factor becomes a step function which is not smooth. Thus, we have 
chosen the minimum n's to achieve the maximum smoothness. The individual contributions 



of the soft {F, \G\} s ° h , twist-2 A{F, \G\}r ' and twist-3 A{F, \G\}r ' are summarized 



in Fig. EJ The soft and the twist-3 terms turn out to give dominant contributions at the low 
and moderate range of Q 2 -values with anomalously large twist-3 contributions in the time- 
like region. Nevertheless, both the corrections exhibit sharp fall-off with increasing Q 2 , such 
that the plot extended beyond Q 2 ~ 50 GeV 2 , will clearly show the twist-2 contributions 
as the being the only dominant ones, both for the space- and time-like domains. Note that 
our factorized hard results are calculated using non-asymptotic collinear twist-2 and twist-3 
DAs (up to NLO in conformal twist), taken from Ball et al. [22[ where these are obtained 
in the context of LCSR. 

The final result for the total space- and time-like form factors (Eq. [TO]) is displayed in 
Fig. El along with the existing experimental pion data. The solid (red) curves correspond 



{F,GUQ 2 ) = 1 



1 + Qs /Q 2 



+ A{F, G}^ ist2 \Q 2 ) + A{F, G}% wist3 \Q 2 ) ; 



(1 + 4s /Q 2 )3/2 




(10) 
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FIG. 3: The total space- and time-like e.m. pion form factors calculated using Eq. (TUJ denoted 
by the solid (red) lines; the soft from factor {i 7 , |G|}| is denoted by the double dashed (black) 
lines; the Ward-identity violating result that does not include the pre-factor modification is also 
displayed. For comparison, the standard asymptotic pQCD result {2\ is displayed. The shaded area 
is roughly our estimated theoretical error beyond Q 2 ~ 5 GeV 2 . The world pion data are taken 



from 



3 61. 



to our central result obtained with the chiral parameter /i n = 1.5 GeV, while the shaded 
area can be regarded as our estimated theoretical error, reliable only beyond the resonance 
region. Due to the overwhelming time-like twist-3 contributions at intermediate energies, 
the modulus of the total time-like form factor \G n \ shows a big enhancement, at least by a 
factor of ~ 2 — 4 compared to F n , although both clearly show asymptotic trends, numer- 
ically approaching the standard pQCD result 2] (beyond ~ 50 — 100 GeV 2 ). The above 
enhancement is largely due to two model independent features: 

(1) The quark and gluon propagators in the hard kernel can become on- shell at non-zero 
q 2 in the time- like but not the space-like region. Thus, in general |G^ ard | should be bigger 
than _p^ ard . In pQCD, this generic feature is captured only if the k? dependence is kept, so 
that the denominator of the hard kernel has terms proportional to q 2 and k% that cancel 

n 

each other in the time-like but not the space-like region (see, e.g., Eqs. 27 and 43 of |I_2j 
before they are Fourier transformed into Eq. [7] of this paper). Without the k? dependence, 
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this enhancement will be missing. 

(2) The twist-3 contributions are more important than the twist-2 ones for intermediate 
range of Q 2 due to two notable reasons: Firstly, the twist-3 terms have the aforementioned 
parametric enhancement arising from fi n which is absent in the twist-2 case. Secondly, the 
finiteness of and the derivative of with respect to x (see, Eq. E]) at the endpoints 
x = 0, 1. These features together with H oc x~ 3 ^ 2 , being oscillatory in the time-like region 
and exponentially decaying in the space-like region (i.e., the time-like parton propagators 
in momentum representation develop poles which are absent in the space-like), account for 
the characteristic relative enhancement of the time-like twist-3 contributions. 

With these two rather robust features and typical treatments of the soft and sub-leading 
twist-3 contributions, it appears that the previously reported discrepancy between the ex- 
perimental data and theoretical predictions jlj] can be ostensibly reconciled. To this end, 
we also present our results without including the pre-factors to demonstrate their effect. As 
revealed from Figs. |2] and EJ without the pre-factors the hard contributions tend to grow 
very rapidly as Q 2 — > and become unreliable, while at the same time beyond Q 2 ~ 5 — 10 
GeV 2 the effect of the pre-factors is hardly discernable. Clearly, then our predictions con- 
vincingly agrees with most of the space- and time-like experimental pion data, including the 
recent CLEO result: Q 2 |GV(13.48 GeV 2 )| = 1.01 ± 0.11(stat)±0.07(syst) GeV 2 and also 



the theoretical prediction M]^ [M 2 ^ = 9.6 GeV 2 ) = 0.94 ± 0.08 GeV 2 [31], fixed from 
branching ratios of J/ip — > titt and J/ip — > e + e~ decays. 

Finally, to comment on the error estimate of our approach, we first look at the error 
band in Fig. |3] While the width of the error band is too narrow to be even noticeable 
in the space-like region, it is anomalously large in the time-like region. Over 90% of this 
error is essentially due to the variation of the chiral parameter fi n between 1.3 — 1.7 GeV 
with increasing contribution to the pion form factor. The remaining difference generously 
over-estimates the other model (parameter) dependences in the DAs from QCDSR, but it 
seems to be a reasonable range of theoretical error when the error of the soft part is also 
included. However, we again stress that the estimate only applies beyond the resonance 
region. Furthermore, several aspects deserve to be noted: The extent of the theoretical 
error from our LO analysis is large enough to completely subsume the systematic errors 



that may arise, e.g., considering NLO effects (in the QCD coupling a s ) [15|, |32j, |33j, sub- 
leading twists (see, e.g., 34j for the twist-4 and twist-6 contributions to the pion form 
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factor in the context of QCDSR), and effects due to higher Fock state corrections which 
are expected to be rather nominal. For example, even without explicit calculations, it is 
easily understandable that the 2-particle twist-4 power corrections are, in fact, very small 
being being proportional to — > 0. Again, the contribution of the 3-particle twist-3 DA, 
being proportional to the "tiny" non-perturbative parameter f 37T « 0.45 x 1CF 2 GeV 2 (to 
be compared with the 2-particle twist-3 DA parameter fi n rs 1.5 GeV), is also strongly 
suppressed. Thus, the 2-particle twist-3 contributions are indeed very special in this regard. 
Moreover, it is estimated that the NLO corrections in the case of the 7T7* — > 7 transition form 
factor amount to only about 5% under specific factorization scheme with the factorization 
scale set to the energy /momentum transfer Q [26j. This is not expected to be very different 
for the pion form factor. To conclude, the unnaturally large twist-3 contribution, especially 
in the time-like region, is certainly non-intuitive and may constitute an important step 
toward understanding the large asymmetry seen in the experimental data, unaccountable 
otherwise. 

The authors thank C.-W Kao, J. Qiu, A. V. Radyushkin, W. Schroers, N. G. Stefanis, 
Z. T. Wei and X.-G. Wu for fruitful discussions. We are especially thankful to H.-N. Li for 
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